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ABSTRACT 

A method for integrated control/structure 
optimization by multilevel decomposition is 
presented. It is shown that several previously 
reported methods were actually partial 
decompositions wherein only the control was 
decomposed into a subsystem design. One of these 
partially decomposed problems was selected as a 
benchmark example for comparison. The present 
paper fully decomposes the system into structural 
and control subsystem designs and produces an 
improved design. Theory, implementation, and 
results for the method are presented and compared 
with the benchmark example. 

INTRODUCTION 

Over the last decade, increasing attention has 
been given to control/structure interaction (CSI) 
problems and the integrated design of a structure 
and its controller. In general, the approaches to 
integrated design can be categorized as either 
sequential or simultaneous [1]. A sequential 
method is one in which one disciplinary design 
iteration precedes the other. The interactions 
between the structure and control are examined in 
an analysis step but are not rigorously included in 
the determination of design changes in the next 
iteration step. A simultaneous design method is 
one in which the control and structure design 
problems are combined into a single design 
problem. Interactions are considered at the outset 
and the effects of the structural and control design 
variables are considered together. References 1 
through 19 give examples of simultaneous 
structure and control design approaches, each of 
them optimization-based. 

Study of these papers shows that there are a 
variety of possible approaches to the problem of 
simultaneous structure and control design even 
within the limits of being based in optimization 
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methods. Often, the objective function is structural 
weight plus some weighted function of controlled 
response. Typical constraints might be on closed- 
loop frequencies or structural deflections. Usually, 
the methods proposed have stated the design 
problem as a single optimization problem wherein 
detailed structural design variables (member sizes) 
and control gains are of equal status as problem 
design variables. The problems on which these 
methods have been exercised have been relatively 
small, consisting of a simple structure with a few 
structural design variables and simple controllers 
with a few gains. While these methods have 
worked well, they could easily require extremely 
large design problems for a large, complex 
structure. To make the problem size more 
manageable, the number of design variables might 
have to be reduced to ineffectual levels. There is 
also an organizational problem associated with 
combining the synthesis of structure and controller 
so directly. A typical aerospace design organization 
is highly segregated along disciplinary lines. 
Obviously, sufficient communication exists 
between disciplinary design groups to effect the 
design of successful aerospace systems. However, 
there are cases of actual flight hardware with 
problems that might well have been caught before 
production had the design process more closely 
coordinated structure and control design functions. 

Multilevel decomposition is an alternative 
approach to large multidisciplinary system design 
that has been proposed [20]. In this approach, a 
large multidisciplinary system is broken down 
along disciplinary and hierarchial lines into 
subsystem designs that are smaller and more easily 
managed than the complete, integrated system. The 
subsystem designs are coordinated at a higher level 
where the influences from each subsystem are 
integrated. At this top level, the design variables 
quantify the Influence of the subsystems on the 
total system, but are not at the level of detailed 
design. Desired changes in these variable are passed 
to the appropriate subsystem design functions 
wherein detailed designs are accomplished treating 
the quantities specified from above as parameters to 
be held constant. An important element of this 
process is the return, from the subsystem designs 
to the upper level, of the values of sensitivity 
derivatives with respect to the parameters held 


constant. These sensitivities are used at the top 
level to construct gradients for use in an 
optimization. The sensitivities may be of any 
quantities that influence the objective and constraint 
functions being used at the top level. If the 
subsystem designs are themselves accomplished 
through optimization, then the sensitivities would 
be the so-called "optimum sensitivity derivatives" 
[21] or, more correctly, "sensitivities of the 
optimum". Thus these sensitivities account for 
how the optimum subsystem design will change as 
the parameters are changed. Further, since the 
controlling influence upon the disciplinary designs 
exists at the top level where the design variables are 
related to each of the contributing subsystems, the 
top level design is a simultaneous one. 

Several of the simultaneous structure and 
control design methodologies cited in this paper [5, 
13, 15-17] have actually been set, knowingly or 
not, in this multilevel format. In these papers, the 
structure was represented completely at the top 
level with no decomposition. The structural design 
variables were either the detailed structural sizes or 
more global structural parameters (such as stiffness 
characteristics) without regard to the detailed 
structural design. All of them, however, also 
incorporated the optimal, steady-state. Linear 
Quadratic Regulator (LQR) as controller. At each 
iteration step, a LQR design is obtained and 
sensitivities of this optimum are computed with 
respect to the top level variables (structural sizes 
and, in some cases, elements of the weighting 
matrices used in the LQR synthesis), and used at 
the top level in constructing gradients of the 
objective and constraint functions for the next 
iteration. Sensitivity equations for the LQR and 
LQG (Linear Quadratic Gaussian for systems 
incorporating Gaussian white noise) control have 
been well developed over the past years [22], 
Sensitivities and even multilevel decompositions 
for purely structural systems under static loading 
have also been developed (23, 24]. However, to 
date, decomposition of both structure and control in 
the controlled structure design problem has not 
been achieved. It is the purpose of this paper to 
report on a simultaneous structure and control 
optimization wherein the total system has been 
decomposed into both structural and control 
subsystems. 

MODEL DESCRIPTION 

The structural model chosen for this study is 
the ACOSS-4 [25], shown in figure 1. The 
tetrahedral truss, mounted on three bipods, is 
intended to model a feed tower for a class of large 
space antennae. Colinear with each bar in the 


bipods, an axial force actuator is assumed to exist. 
All geometric and material properties are assumed 
to be in unspecified but consistent units. Table I 
gives the nodal coordinates, and Table II gives the 
nodal connectivities and the initial cross-sectional 
areas of the truss elements. The elastic modulus 
of each member is 1.0 and the mass density is 
0.001. A concentrated mass of 2.0 is assumed to 
be at each of the four comers. Further, the mass 
"units" are assumed to be such that the total mass 
of the structural elements is negligible in 
comparison to the total of the concentrated masses. 
Thus the structural mass does not enter into the 
dynamic characteristics of the model. This model 
is a good compromise between simplicity and 
complexity, and is well represented in the open 
literature. The controller is chosen to be the 
optimal steady-state LQR. Its use will also allow 
comparisons to be made with published results 
[17]. 


PROBLEM 

Equations of Motion and LQR Control 

The matrix equations of undamped, forced 
motion for the Finite element model of the structure 
take the form, 


Mw + Kw = bu (1) 

where M and K are the mass and stiffness matrices, 
and w is the vector of x, y, and z displacements of 
each of the comer nodes. The matrix b is the input 
matrix defining the geometry of the applied load 
distribution. In this case it contains the direction 
cosines of the six actuators. The input vector u, 
then, is composed of the forces exerted by each of 
the actuators. The undamped natural modes of the 
structure are obtained by solving the eigenvalue 
problem, 

K - eo? M ] i = 0 (2) 

Assuming linearly independent eigenvectors, the 
damped equations of motion can then be cast in 
modal form, 


M x + 2D. £x + Kx = b u , (3) 

where 

M = d> T M<r>=i 

k =<t> T K<i> = n 2 

b = <t> T b 
and 

w=<bx 
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Above, Q is a diagonal matrix of modal damping 
factors, <J> is the modal matrix containing the 
eigenvectors obtained in the solution of equation 2, 
is a diagonal matrix of squared structural 
frequencies from the solution of equation 2, x is 
the vector of generalized coordinates, and I is an 
identity matrix. 

The equations of motion may then be put in 
state-space form, 

X = AX + Bu (4) 

where 

x=U} 

0 I 

■si 2 -2SIC, 


A linear state feedback control law would be of 
the form, 

u = *GX (5) 

The LQR control law is that which minimizes the 
cost function, 

J=|{ X T QX + u T Ru}dl (6) 

o 

subject to the constraint of equation 4. The 
matrices Q and R are weighting matrices. The 
gain matrix is found from, 

G = R''b T P (7) 

where the symmetric matrix P is the minimum 
solution to the steady-state, algebraic matrix 
Riccati equation, 

A T P + PA-PBr'b T P + Q = 0 (8) 

Then, the equations of motion of the closed-loop 
system may be written, 

X = A cl X (9) 

where 

= A - BG 

Benchmark Example 

The design problem to be considered, taken 
from reference 17, is the minimization of the 
Frobenius norm of the feedback gain matrix, 
subject to lower bound (lb) constraints on the first 
two closed-loop (cl) frequencies and the damping 
parameter of the first closed-loop mode. Side 



constraints are imposed on the cross-sectional areas 
of the truss members, which are the design 
variables, y. Also, the open- loop damping factors, 
are assumed to be zero. The problem reads: 


subject to: 


min Fr = Tr (G R g G) 


0) >CiK 
*\ d l 


CO >00, 


- ^cl 
1 1 

y lb £y<y ub for each y 


(10a) 

(10b) 

(10c) 

(10d) 

(10c) 


where, for eigenvalue 


X= o±jto 

the damping parameter is defined. 


~2 2 

V o +CD 

The gain weighting matrix, R 0 , is taken to be an 
identity matrix and is distinct from the control 
weighting matrix in equation 6. The weighting 
matrices Q and R in equation 6 are also identity 
matrices. In reference 17, the lower bounds in 
constraint equations lOb-d are 1.341 r/s, 1.6 r/s and 
0.15, respectively. The lower bound on each of the 
design variables was 10 while two different upper 
bounds (ub) were considered. One was 1000, and is 
referred to as ’Design IT and the other was 2000, 
referred to as Design UT. 

A description of the method used in reference 17 
will be given partly because similar ideas are 
important for the present formulation. It will also 
facilitate demonstrating how the approach of 
reference 17 (as well as references 5, 13, 15, and 
16) was actually a partial decomposition. Some 
results from reference 17 will be shown later for 
comparison with the present results. Gradients of 
the Frobenius gain norm and the constraints are 
found from. 
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guarantees an asvmtotically stable closed-loop 
system, the coefficient matrices (which are A ci and 

A c1 t ) of dP/dy are asymtoticaily stable. Thus an 
unique solution for dP/dy is guaranteed [22]. 

As mentioned earlier, this organization of the 
design problem is a partial multilevel 
decomposition, as shown in figure 2, wherein only 
the control design is decomposed into an optimal 
subsystem design. The determination of the 
detailed structural design variables is handled at the 
top level and there is no structural decomposition. 
The structural design variables serve as constant 
parameters for the control design and sensitivities 
of the optimal control solution, and resulting 
behavior variables (closed-loop eigenvalues), are 
computed. As described in reference 21, the 
optimum sensitivities are obtained through 
differentiation of the conditions of optimality 
satisfied by the solution of the optimization 
problem. The Riccati equation (equation 8) fills 
the role of optimality conditions for the optimal 
steady-state LQR problem. 

Present Formulation - General 


* 

where each diagonal element of the derivative of the 
£5 2 matrix is, 

3co j 2 TfiK n :3M ) 

5 y = V dy j dy jv ^ 

In the interests of brevity, the expressions for the 
derivatives in equation 16 of the eigenvectors will 
not be reproduced here as they arc given in reference 
17. 


Finally, the key element is the derivative of the 
feedback gain matrix, 

rp- = R P +■ R"'b T -IM- (18) 

dy dy dy 

The required derivative of the Riccati solution, P, 

is obtained by differentiating the Riccati equation to 

yield, 

f \ ( N'T 


A - Br 'b T P 


+ P-l^+^r-P P 

dy dy 


A -BR 'b T P 


3^-R V+BR ' 1 -^- 
dy dy 


This equation is a Lyapunov equation with 3P/3y 
as the unknown. Since the LQR control solution 


A generalized decomposition would appear as 
shown in figure 3 wherein the objective function is 
minimized with respect to higher level design 
variables, denoted by z’s, that characterize 
subsystems I and II. In each iteration, the 
subsystem optimizations are performed with 
equality constraints enforcing the invariance of the 
parameters, z, specified as design variables at the 
higher level. After each subsystem optimization, 
sensitivities of the optimum solutions are 
computed, as are sensitivities of appropriate 
behavior variables, denoted in figure 3 by q's, that 
may be of importance to the top level design. The 
behavior variables might be used in the top level 
objective or constraint functions. In the 
development that follows, optimum sensitivities 
are indicated by superscript asterisks. 

The general top level optimization problem is 
given by, 

min F(z, q) (20a) 

subject to: 

g.(z, q) < 0 (20b) 

^Ib - Z — Z ub (20c) 

The optimization problem for subsystem I would 
be, 

min F, (y, . q, (y, ), q n (y n ) ) (21a) 
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subject to: 

g[ ,(y r qi (y, ). q n (y n > ) s 0 ( 21b > 

Zi=z I (y I .q,(yi).q u (y n )) ( 21c ) 

y V Sy|Sy 'uS (21<1) 

tt> ub 

Note that dependencies upon behavior variables, q n , 
from subsystem II are included in the functional 
statements above. This is because there may be 
strong physical coupling between the two 
subsystems. If so, then a certain amount of 
iteration between the subsystem designs may be 
necessary at the lower level. The first constraint 
represents those constraint conditions applicable 
directly to the subsystem design. The second 
constraint represents the requirement that the top 
level design variables applicable to subsystem I be 
the parameters held constant, or matched, during the 
subsystem optimization. The statement of the 
optimization problem for subsystem II is similar. 


Gradients of the objective and constraints for 
the top level optimization are, 


5a* 

dF = ill + V 1L _L 
dz 3z ^ dq ; 5z 

j J 


( 22 ) 


and 

(a) 

dz dz 2mJ 3q 5z 
j J 

The 5/5z derivatives are pseudo-partial derivatives 
in that while they act as partial derivatives in 
equations 22 and 23, they are constructed from 
sensitivities from the sublevels. The superscript 
asterisks indicate that optimum sensitivities from 
the subsystem designs are used, 



The above equations are quite general and simplified 
in that no distinction is made between q’s, y's, and 
z's for one subsystem and those for another. 
However, to do so would unnecessarily complicate 
the presentation of the equations and would not 
increase their generality. Application of these 
equations to any particular problem would involve 
careful attention to actual functional dependencies 
among the problem variables and parameters that 
would be unique to each application. The optimum 
sensitivities in the above equations come from 
differentiation of the necessary conditions of 
optimality, derived in reference 21 in a form that is 
typical of structural optimization problems. The 
essential information for computing the 
sensitivities of the optimal steady-state LQR 
problem is given briefly in this paper and in more 
detail in reference 22. Accordingly, further 


analytical description of these sensitivities will not 
be given here and the interested reader is encouraged 
to consult references 21 and 22. 


Present Formulation - Detailed 


The distinguishing element of the present 
formulation will be the decomposition of the 
structural subsystem to a sublevel equal to that of 
the control subsystem. Since the problem to be 
solved is the same as the benchmark example, the 
top level objective and constraints are similar to 
those in equations lOa-d, except that different top 
level design variables will be used. Also, the 
constraints are formulated as reciprocal Taylor 
series approximations. Since the benchmark 
example only uses structural design variables at the 
top level, the new top level design variables will 
only be structural as well. The structure is 
represented in the control design by its modal 
characteristics, so the top level design variables 
will be a subset of the squared natural frequencies of 
the open-loop structure. Specifically, the squared 
frequencies of modes 1, 2, 4, 6, 8, 10, and 12 will 
be used. In the following discussions, these will 
be denoted by a> p 2 *s, corresponding to the z r s in 
the general formulation given previously. The 
remaining squared frequencies will be denoted by 
co N 2 's and will be treated as behavior variables, 
q s 's, from the structural subsystem. 


The statement of the top level optimization 
problem is, then, 

(25a) 


min Fr = Tr (G R fi G) 
2 


subject to: 


- V ** 
8 > = g 0 + 2-TT 
‘ > 


7 

( \ 
2 2 

V\ 


\ ) 


<;o (25b) 




(25c) 


where, for a> c | , 


8 s 


- < 0 , 

cl, d. 


CD . 


The subscript 0's denote starting values. The 
expressions for (0 c i 2 and £ c j, are similar. Move 

limits of ±p of the current starting values at each 
step are imposed on the design variables as side 
constraints. 
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The structural subsystem design is posed as an 
optimization problem to minimize weight with 
respect to the cross-sectional areas of the twelve 
truss members subject to the parameter matching 
constraints (i.e. equation 21c) and side constraints 
on the design variables. The authors did not wish 
to attempt enforcing strict equality for the 
parameter matching constraints, so the 0 ) p 2 s 
desired from a top level design iteration were 
considered matched if the corresponding eigenvalues 
from the structural design were within a small 
tolerance, ±e. Thus the problem reads. 


nyn W(y) 


(26a) 


subject to, for each desired parameter, : 

( , xV 

«\> 2 (y) 


i - 


1 


% 


doB 


1 < 0 


and for each variable, y, 

10* y £ 2000 


(26b) 


(26c) 


There are no additional constraints (the gfs of 
equation 21b) such as stress or deflection limits 
since there are none imposed in the benchmark 
example. 


The LQR control subsystem design is the same 
as described by equations 5 through 8. The real and 
imaginary parts of the closed-loop eigenvalues are 
behavior variables, q^., from the control subsystem. 
The expressions for the sensitivities are similar to 
equations 1 1 through 19 except that y's are replaced 
by G) p 2 's for parameter sensitivities and by co N 2 ’s 
for structural behavior variable sensitivities. 
Sensitivities of B are zero since it is assumed not a 
function of frequencies. 


Once the subsystem optimizations arc 
complete, optimum sensitivities with respect to 
parameters and behavior variables can be computed 
and the top level gradients can be constructed. Note 
that the physical coupling between the structure and 
control designs is one way, namely, from structure 
to control. Thus, none of the control behavior 
variables influences the structure design, whereas 
structural behavior variables (squared frequencies) do 
influence the control design. The first step is to 
calculate the pseudo- partial derivatives of behavior 
variables in equation 24, For the structural 
behavior variables, qj = co^ 2 , equation 24 becomes. 



The derivatives of the squared frequencies with 
respect to the structural design variables, y, are 
available from the implementation of the structural 
optimization step and are computed from equation 
17. While they are not optimum sensitivities 
because they are not in any way computed from 
differentiation of the optimality conditions, they are 
the eigenvalue sensitivities for the structure at the 
optimum solution. For the control design behavior 
variables, equation 24 becomes (for the closed-loop 
eigenvalue imaginary parts only). 


5<d^ 5g) c1 9co d 

j 


fcflp d<a p 


hi 


„ 3< \ V 

*4, C8> 


Note that the last summation in equation 28 is 
simply the expression in equation 27. The 
expression for the eigenvalue real parts is similar. 


The gradients for the top level objective and 
constraint functions can now be obtained, with 
reference to equations 27 and 28, as 


dFr 


( T * 

aG R_G + G T R, 


A 


V 




2 G' 


0,2 


J 


f 


\ 


SG R..G + G T R iC 


. 2 G 

3t0 N 
' j 


G -* 2 


5{u n 

* j 

Scop 2 


(29) 


dgj dg. 


dc£ do^ 


(<) 

UJ 


where, For co cll 


(“d ) 3 I(“, ) Sco* 


1 l 




^ H 2 


(30) 


Implementation 

The top level and both sublevel optimizations 
were written as separate codes that communicated 
through data files. The ADS [26] optimization 
code was used at the top level and structural 
subsystem while the optimum LQR solution was 
obtained using the ORACLS optimal control 
design software package [27]. The Sequential 
Quadratic Programming strategy was used at the 
top level with the Modified Method of Feasible 
Directions optimizer and the one dimensional 
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search option wherein the minimum of the 
constrained function is found by first finding the 
bounds and then using polynomial interpolation. 
The same options were used in the structural 
subsystem optimization. The move limit, |i, on 
the top level design variables was set initially to 
5%. This was reduced as the final solution was 
approached to overcome oscillation in the 
constraints. The tolerance, e, in the parameter 
matching constraints was set to 1%. 


The structural subsystem optimizations tended 
to experience difficulty satisfying the parameter 
matching constraints (equation 26b) since any top 
level design variable change of more than 1% 
resulted in an infeasible starting point at the 
structure sublevel. Thus, an intermediate step was 
introduced into the procedure by which the 
structural design variables were adjusted until the 
1% parameter matching constraints were met. 
These values of the structural design variables then 
served as the starting point for the structural 
subsystem optimization. The adjustments were 
automated by forming an optimization problem 
(again, using ADS) wherein a cumulative function 
of the parameter matching constraints is minimized 
with respect to the structural design variables until 
the function becomes negative, thus ensuring that 
the constraints were satisfied. A cumulative 
function with this desirable property is the so-called 
KS function [28]. For a set of NC constraints, g v 

the function is , 

NC 

= hi^e (31) 

i- l 

where p is a weighting parameter selected so as not 
to produce an over- or underflow condition when 
computing exponentials of large positive or 
negative numbers. This function has the property 
that it is always greater than or equal to the 
maximum constraint within a definite band. 


f In (NC) (32) 


The algorithm proceeds as follows. The 
matching tolerance is set to 4% (one percentage 
point less than the top level move limit) and the 
first matching constraint is put into the cumulative 
function which is reduced by the optimizer until it 
becomes negative. Then the first matching con- 
straint serves additionally as an actual constraint 
while the second matching constraint is added to the 
cumulative function which, again, is reduced by the 
optimizer until it becomes negative. This 
continues until all of the matching constraints 
(seven of them, in this case) satisfy the 4% 
tolerance. The tolerance is then reduced to 3% and 
the process starts again with the first matching 


constraint. This continues until all of the 
matching constraints satisfy the 1 % tolerance. 

NUMERICAL RESULTS AND 
DISCUSSION 

Figures 4 and 5 show the line-of-sight (LOS) 
error responses (total lateral deflection of node 1 
away from the z axis) of the open- and closed-loop 
initial structure to a unit displacement of node 2 in 
the x direction imposed as an initial condition. The 
control manages to suppress the large excursions, 
but reduces the average deflections little. Also, the 
damping appears light to moderate. In contrast, 
figure 6 shows the response of the present design. 
The overall response is drastically reduced as are the 
larger excursions, and the system is more highly 
damped. Note the expanded scale of figure 6. 

Structural member sizes, and open- and closed- 
loop eigenvalues for the present design are shown 
in Tables III through V along with corresponding 
results for Design IF and Design nr from reference 
17 and the initial design. It can be seen in Table 
IV and figure 7 that the present design is an 
improvement over Designs II and in’ in not only 
the Frobenius gain norm, which was the objective 
function in all cases, but also in the weight of the 
structure, which was the objective function in the 
present method's structural subsystem design. 
Whereas, in the benchmark examples, the member 
sizes were allowed to range widely so as to satisfy 
the problem constraints as easily as possible, they 
were rather more constrained in the present example 
to produce the desired open-loop frequencies alway s 
at a minimum weight. This demand is undoubtedly 
responsible for the reduced weight of the present 
design. Since the structural mass is assumed 
negligible in the dynamics, it is consistent with the 
reduced weight that both the open- and closed-loop 
frequencies in general show the present design to be 
a more flexible one. 

Several attempts were made to reduce the degree 
to which the physical coupling between the 
subsystems was accommodated. One approach was 
to pass the modal information resulting from the 
parameter matching program to the control design. 
The other was to pass the desired frequencies from 
the top level and the mode shapes from the 
previous structural design. Both approaches failed 
to converge to a feasible design. It is believed that 
the sensitivities from the subsystem designs used at 
the top level were inconsistent with one another 
resulting in erroneous top level gradients. 
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CONCLUSIONS 

A method for integrated control/structure 
optimum design by multilevel decomposition has 
been presented. It was applied to a design problem 
found in the literature and the results were 
compared. The method was found to produce lower 
values of the objective function, which was the 
Frobenius norm of the state feedback gain matrix, 
as well as the structural weight. It is concluded 
that using structural weight as the decomposed 
structural subsystem design objective has a 
moderating influence on the weight of the design. 
It was also found that the physical coupling 
between the subsystems mus t be con sidere d in the 
design procedure. The structurc-to-control coupling 
that is present in this example will occur whenever 
the control synthesis is model-based. 
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Element 

Nodal 

Connectivities 

Cross Sectional 
Area 

1 

1 - 2 

1000 

2 

2 - 3 

1000 

3 

2 - 4 

1000 

4 

3 - 4 

1000 

5 

1 - 3 

100 

6 

1 - 4 

100 

7 

2 - S 

100 

8 

2 - 6 

100 

9 

3 - 7 

100 

10 

3 - 8 

100 

11 

4 - 9 

100 

12 

4 - 10 

100 


Table III ACOSS-4 Member Cross- 
Sectional Areas 


Mew her 
No. 

Initial 

Design 

'Detiga n 

(Ref. 17) 

Design ni* 
(Ref. 17) 

r reseat 
Desiga 

1 

10M 

990.07 

175.20 

1040.10 

2 

10W 

990.07 

1991-30 

1014.40 

3 

10W 

990.07 

1351.59 

1054.02 

4 

1000 

990.00 

1034.27 

1023.21 

5 

100 

990.58 

1902.14 

171.0 

4 

100 

990.10 

1905.75 

242.75 

7 

100 

921.04 

102.38 

140.11 

3 

100 

902.75 

159.91 

153.71 

9 

100 

30.93 

29.20 

90.99 

19 

190 

194.24 

152.80 

217.97 

11 

100 

30.09 

29.19 

29.00 

12 

100 

07.51 

407.09 

201.02 

Froheaiu* 

Norm 

17.175 

21.909 

18.01 

17.945 

Weight 

43.49 

00.080 

I7i34 

47,843 


Table I ACOSS-4 Nodal Coordinates 


1 

X 

V 

z 


0.0 

0.0 

10.165 


-5,0 

-2.887 

2.0 


5.0 

-2.887 

2.0 


■m 

5.7735 

2.0 


-6.0 

-1.1547 

0.0 


-4.0 

-4.6188 

0.0 


4.0 

-4.6188 

0.0 

* 

6.0 

-1.1547 

0.0 


-2.0 

5.7735 

0.0 


2.0 

5.7735 

0.0 


Table IV ACOSS-4 Structural Eigenvalues 
(Squared Frequencies) 


Mode 

No. 

Initial 

Design 

'Design IT 
(R«r. 17) 

'Design 111 1 
(Ref. 17) 

Present 

Design 


1.80 

2.08 

If 

mm 

2 

2.77 

3.13 

1 


3 

8.35 

5.57 




8.74 

15.91 

1531 

WEEM 

5 

11.54 

33.49 

17.00 

17.72 

0 

17.07 

3931 

34.95 

1838 

mm 

21.73 

85.79 

42.95 

38*83 

B 

22.01 

94.51 

9037 

39.49 

B 

72.92 

101.50 

112.70 

81.19 

kJ§ 

85.57 

17730 

17230 

9439 

Gfl 

105.77 

202.70 

214.00 

111.9* 

H 

100.54 

23130 

310.90 

174.75 
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Table V ACOSS-4 Closed-Loop Eigenvalues 



'Design II' 
(Ref. 17) 

'Design III' 
(Ref. 17) 

Present 

Design 

- 0.073 ± 1.342j 

- 0.246 ± 1.435j 

- 0.252± 1.369j 

- 0.207 ± U42j 

- 0.109 ± 1.663j 

• 0.256 ± 1.762j 

- 0.251± 1.670j 

- 0.126 ± 2.1 78j 

- 0.213 ± 2.885j 

- 0.274 ± 2.346j 

- 0.185 ± 2.1 13j 

- 0^32 ± 2J72j 

- 0.237 ± 2.950j 

- 0.266 ± 3.980j 

• 0.267 ± 3.905j 

- 0220 ± 3.720j 

- 0.285 ± 3.388j 

- 0.275 ± 5.781 j 

- 0295 ± 4.185j 

- 0^62 ± 4.202j 

- 0.363 ± 4.191j 

- 0.231 ± 6.265j 

- 0.344 ± 5.903j 

- 0J37 ± 4.298j 

- 0.355 ± 4.649j 

- 0.236 ± 9.259j 

- 0.265 ± 6.548j 

- 0J54 ± 6.222J 

- 0.344 ± 4.744j 

- 0.261 ± 9.718j 

- 0209 ± 9.812j 

- 0J02 ± 6.239J 

- 0.292 ± 8.535j 

- 0.180 ±10.075j 

- 0.290 ± 10.6l3j 

- 0-300 ± 9.005j 

- 0.276 ± 9.247j 

0.217 ± 13.3 14j 

- 0.292 ± 10.004j 

- 0.284 ± 9.737j 

- 0.214 ± 10.283j 

0.339 ± 14,234j 

- 0.231 ± 10.777j 

- 0.222 ± 10.58j 

- 0.083 ± 12.905j 

0.284 ± 15.200j 

- 0.092 ± 13. 323 j 

- 0.091 ± 13.219j 



Figure 1 ACOSS-4 Structural Model 
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Figure 3 Generalized Decomposition 
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Figure 5 Closed-Loop Response of 
Initial Structure 
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Figure 6 Closed-Loop Response of 
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